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Numerical  Model  of  Dispersion  of  Atmosphere  Pollution  Within  Urban 
Agglomeration  .  (M.  Uliasz) 


Introduction. 


One  of  the  most  frequent  of  air  pollutants  emitted 


in  urban  environments  is  sulfur  dioxide.  Oxidation  processes  in  the  atmos¬ 
phere  are  the  cause  of  sulfur  dioxide  turning  into  sulfuric  acid.  As  a  result 
of  chemical  reactions  with  1  mg  of  SOg,  1.  531  mg  of  H2SO4  occurs.  The 
acceptable  concentrations  for  safety  zones  are  about  3  times  stronger  for 
sulfuric  acid  than  for  sulfur  dioxide,  or  that  same  quantity  of  sulfur  in  the 
form  of  H2S04  in  free  air  is  about  a  4.  5  times  greater  threat  to  the  environ¬ 
ment  than  in  the  form  of  sulfur  dioxide. 


The  paper  presents  a  numerical  model  of  the  dispersion  and  conversion 
of  sulfur  dioxide  emitted  from  an  urban  area.  Using  a  two-dimensional  diffu¬ 
sion  model  makes  it  possible  to  find,  besides  the  SO2  decomposition  concen¬ 
tration,  the  decomposition  concentration  of  sulfuric  acid  also  as  one  of  the 
end  products  of  the  oxidation  of  sulfur  dioxide. 


Diffusion  Model. 


given  by  the  equation 


The  model  considered  in  pollution  dispersion  is 


*  .LgU  +  f 
"  TS  it  *  is  r 


1. 


with  boundary  conditions 
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S4  s  (x,  z)  is  the  pollution  concentration,  fj  is  the  volumetric  intensity  of 

G-f?ouMj> 

the  pollution  source,  v^  is  the  pollution  absorption  rate  by  the-dant,  Qj  is  the 
pollution  emission  from  a  unit  of  surface,  whereby  when  i-  1  the  quantity 
reaches  SOg,  and  when  i«  2  the  quantity  reaches  H2SO4,  usu(z)  the  velocity 
of  air,  K=K(z)  the  vertical  turbulence  coefficient,  H  is  the  height  of  the  mix¬ 
ing  layer . 


Fig.  1.  Field  and  boundary  conditions  of  the  model. 


The  field  of  action  of  the  model  is  presented  in  Figure  1,  assuming 


that  from  the  surface  of  the  city  sulfur  dioxide  is  emitted  and  with  that,  it  is 
assumed  that  pollution  does  not  take  place  on  the  windward  side  of  the  source 
of  emission. 

Where  zsH,  a  second  variant  of  the  boundary  condition  was  tested  also: 


as 
a  t 


«  o 


With  the  proposed  height  of  the  mixing  layer  H=  1, 000  m,  the  results  from 
the  model,  with  different  boundary  condition  variants  (4)  and  (5),  were  prac¬ 
tically  identical. 


The  model  takes  into  consideration  the  pollution  absorption  by  the 
ground.  The  absorption  speed  is  understood  as  the  relationship  of  the  verti¬ 
cal  turbulent  pollution  flow  to  its  concentration  at  the  surface  of  the  gro  und. 
It  is  assumed  that  sulfuric  acid  is  absorbed  by  the  ground  identically  like 
sulfur  dioxide,  consequently  Vj=V2.  In  open  space  it  is  granted  that  the 
absorption  rate  is  equal  to  1  cm  /s  / 3 / ,  while  for  the  city  it  is  convention¬ 
ally  3  times  greater. 

The  two-dimensional  diffusion  model  defines  the  decomposition  of 
pollution  concentration  emitted  by  the  surface  source  infinitely  extending  in 
the  direction  of  the  y  axis.  In  the  case  of  the  town--a  source  having  finite 


measurements- -the  model  describes  the  concentration  only  in  the  zone  in 
which  the  concentrations  are  not  dependent  on  the  yQ  dimension  of  the  emis¬ 
sion  source  (the  area  is  marked  off  in  Fig.  1).  This  zone  is  designated  by 
the  inequality 

.  (yjz  -  c  sfT~~Z0)  <  y  iyjl  c  n/*  -  *0)  '  (#) 

i  .  . 

where  f  *  3%/ 2  K'Ju,  Ky  -  the  turbulence  coefficient  in  the  direction 
of  the  y  axis. 

The  length  of  the  validity  zone  of  the  model  in  the  direction  of  the  wind 
is  usually  about  4  to  5  times  greater  than  the  width  of  the  emission  source  of 
yQ.  The  model  gives  the  highest  values  of  pollution  concentration  from  those 
that  take  place.  Beyond  the  limits  of  the  field  defined  by  condition  (6)  the 
results  of  the  model  will  be  exaggerated. 

In  the  model  considered  it  is  assumed  that  the  only  product  of  the 
sulfur  dioxide  oxidation  is  sulfuric  acid.  The  SO2  loss  in  equation  (1)  is 
shown  by  fj,  while  the  appearance  of  H2SO4  is  taken  care  of  by  fg,  whereby 
we  get:  fg®  —  1.  531fj.  Most  often  it  is  agreed  that  the  loss  of  sulfur  doxide 
is  proportional  to  its  concentration,  so  called  f^s— kS^.  In  this  case  the  time 
of  half  loss  of  pollution  is  Ti  / 2  =  ln2/k.  In  the  above  supposition  the  calcula¬ 
tions  for  small  wind  velocities  give  relatively  high  values  of  H2SO4  concen¬ 
trations  at-  great  distances  from  the  source  of  emission,  even  higher  than 


from  SC>2  concentrations  A 7 •  Under  actual  conditions  such  situations  are 


not  observed.  In  accepting  the  hypothesis  that  the  conversion  rate  of  SOg 
into  H2SO4  depends  on  the  sulfur  dioxide  concentration,  the  model  uses  the 
following  form  of  conversion  function 


/,  «  +  k2s ]  *  *3  Sj) 


In  this  supposition  the  time  of  half  loss  of  SO2  is  a  function  of  its  concentra¬ 
tion: 


rl/2 


In  2 

*,  ♦  *2Si  +  Ml 


The  conversion  factors  kj,  k2«  kg  were  selected  in  order  to  retain  proportions 
between  the  concentration  of  H2SO4  and  SOg,  occurring  in  actuality. 


In  accepting  that  •  4  •  10*  S~l,k2  *  2  •  10~4  S~*  m3 mg"1  Hj  » 0,1  S'*  m* 

3 

for  the  concentration  0.  04  mg  /m  ,  ^1/ 2  =  ,  and  for  the  concen¬ 

tration  Sj*  0.  004  mg  /m^^y2*^^rs* 

The  numerical  solutions  of  the  diffusion  equation  make  it  possible  to 
take  into  consideration  arbitrary  distributions  of  meteorological  parameters. 
The  wind  velocity  and  turbulence  factor  are  accepted  in  the  model  as  functions 
of  height  over  the  surface  of  the  earth  coarseness  of  the  base  zQ  and  geo- 
strophic  wind  velocity  G: 
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The  dynamic  velocity  u#  can  be  expressed  by  the  Lettau  formula  [A] : 


u.  -  G 


0,16 


Ig  Ro  -  13 


where  Ros=G/(z0l)  is  Rossby's  number,  t  is  the  Coriolis  parameter,  u  is  the 
Karman  constant. 


Fig.  2.  Distribution  of  wind  velocity  and  turbulence  coefficient  in 
the  city  and  beyond  the  city  for  G=6  m/s. 


For  the  urban  area  the  parameter  of  courseness  is  zQ  *1.  5  m  ,  and 
for  the  non-urban  area  z  s0,  03  m.  A  fall  in  wind  velocity  and  an  increase  in 
turbulence  factor  is  attendant  to  an  increase  in  the  coarseness  of  the  ground 
zQ  (Fig.  2.). 
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From  the  equation  for  continuity  for  the  air 


ii  ix 


we  conclude  that  at  the  city  limits,  with  a  change  in  profile  of  wind  velocity  u, 
vertical  movements  appear  which  surround  the  air  with  the  velocity  w.  In  the 
model  these  movements  were  omitted  since  the  value  of  the  vertical  air 
velocity  w  does  not  exceed  1%  of  the  wind  velocity  u. 


Equation  (1)  is  solved  numerically  by  a  network  method  using  the  fol¬ 
lowing  differential  scheme  [6] : 

s/+i.*  ~s/k  Kk+m(sj*i*+i  ~ ~ ~5/*i.k-0  A 

Uk— — — —  «  -I—'  '  . . 

A  Jr  A  *2  1  •. 


when  k=2,  .  .  .  ,  n  -  1,  jsl,  .  .  .  ,  m  -  1. 

Boundary  conditions  (2),  (3),  and  (4)  are  approximated  by  the  following  dif¬ 
ferential  equations: 


\  ,  -  0 


A  *  I.  ...  ii 
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A  +  '  =  2 . 
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To  solve  the  differential  scheme  (12)  with  boundary  conditions  (13)  the 
"progonki"  method  C$]  is  used  which  permits  passing  from  quadratic  equa¬ 
tions  (14)  to  a  linear  equation  system  (15). 


7. 


Ak.i  +  \  Sk  1 ,/  ♦  1  Ck.j*  I  Sk.)+  t  + 

*  *  2,  ...  n  1;  /  =  1,  .  ..  m  -\ 

For  the  given  j  the  solution  sought  is  in  the  form 

Sk  ~  Xk* 1  Sk+ 1  +  Yk+l  k  ~  n ■'* .  1 

where^j  *  p\>  Yi  s  R\ 

Bk  Ak  Yk  +  Fk 

xk+\  *  c  -  x  a  1  *  c"-T"a  k  *  2 . n~I 

and  =  (*2  +  /‘j  yj/d  -  P2  Xn) 

The  "progonki"  method  is;stabile  if 

Ak  >  0,  Bk  >  0  i  Ak  *  Bk  <  Ck 

Analysis  of  Results.  The  model  was  tested  for  3  variants  of 

geostrophic  wind  velocity  3.  6  and  10  m/s.  Uniform  emission  of  sulfur  diox- 
ide  from  the  urban  area  was  assumed,  equal  to  100  t/knr  /per  yr. 

Distributions  of  sulfur  dioxide  concentrations  for  different  variants  of 
geostrophic  wind  velocity  are  presented  in  Fig.  3a,  b,  c.  An  increase  in  wind 
velocity  causes  a  definite  fall  in  the  concentration  of  pollution.  With  an  in¬ 
crease  in  wind  velocity  individual  isolines  of  concentration  enclose  smaller 
and  smaller  areas.  Maximum  concentrations  of  sulfur  dioxide  occur  at  the 
earth's  surface  at  the  leeward  outskirts  of  the  city  (z=0,  x=10  km.  ). 
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Fig.  4  Isolines  of  H2SO4  ^concentrations  [nlm*J 
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Relationship  of  quantity  of  sulfur  in  the  form  of 
H2SO4  to  summary  quantity  of  sulfur  when  0  wf# 


The  effect  of  the  increase  in  wind  velocity  on  the  fall  of  sulfuric  acid 
concentration  is  even  stronger  than  in  the  case  of  SOg  (Fig.  4a,  b,  c). 

The  greatest  H2SO4  concentrations  occur  at  a  certain  height  over  the 
earth  and  are  clearly  displaced  in  the  direction  of  the  wind,  in  relation  to  the 
greatest  concentrations  of  sulfur  dioxide.  An  evident  fall  in  the  concentra¬ 
tions  of  both  pollutants  at  the  earth's  surface  (beyond  the  emission  area  in 
the  case  of  SO2)  is  connected  with  the  absorption  of  the  pollutants  by  the  ground. 


11. 


The  ratio  of  the  quantity  of  sulfur  in  the  form  os  H2S04  to  the  sum¬ 
mary  quantity  of  sulfur  when  zsO  retains,  beyond  the  outskirts  of  the  city, 
values  approximately  constant  for  the  given  wind  velocity  (Table).  Benarie, 
Nonat,  et  al.  A7  show  that  the  values  of  this  relationship  are  weakly  depend¬ 
ent  on  the  distance  from  the  emission  source  and  fluctuate  from  5  to  23%. 
According  to  tests  made  by  Hryniewicz  [2]  for  Warsaw,  the  percent  of  sulfur 
in  a  sulfur  aerosol,  calculated  in  relation  to  the  summary  quantity  of  sulfur 
in  gas  and  aerosol,  amounts  to  7  to  32%  depending  on  the  season  of  the  year. 
The  results  of  the  model  presented  in  the  paper  are,  consequently,  concurrent 
with  the  experimental  data. 


Typescript  published  3/29/77. 
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